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¢ Top Ranked large scale space mission in 2010 New Worlds, New Horizons 
Decadal Survey for Astronomy and Astrophysics 


— Measures Dark Energy, Exoplanet Microlensing, and the near InfraRed Sky 


¢ Includes a 2.4 m existing telescope donated from elsewhere in Federal 
Government 


¢ Includes two baseline instruments supported by Instrument Carrier 
— Wide Field Instrument (WFI) with 2 channels 


— IR imaging with 3x6 array of H4RG detectors for a FOV about 100x 
Hubble’s WFC3 Instrument 


— Integrating Field Channel using a slicer and spectrograph to provide 
individual spectra of each slice 


— CoronaGraph Instrument (CGI) 


— Imaging and spectroscopic modes to image exoplanets and debris discs 
around nearby stars 
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¢ Orbits Earth-Sun Lagrange Point 2 
¢ Spacecraft bus (SC Bus) provides power, attitude 
control, comm., and other spacecraft functions 
— 7 modular, on-orbit serviceable avionics bays” 
¢ SC Bus Top Deck supports 
—I/Instrument Carrier via 3 bipods 
—Solar Array Sun Shield (SASS) to provide 
stable thermal environment 
¢ Instrument Carrier (IC) supports 
— Telescope 
— Wide Field Instrument (WFI) em 2 
— CoronaGraph Instrument (CGI) “he NCGI ace | 
— Instruments are serviceable on orbit 0 — “| Wey’ 
¢ Outer Barrel Assembly (OBA) mitigates stray 
light for telescope. Supported by bipods to SC 
¢ Joint mission by GSFC (BUS, WFI), JPL (CGI), and 
Industry (Telescope and WFI) 
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¢ Survey type missions fall into generally one of two categories 


— (1) Dedicated instrument with large field of view/low resolution searches for source 
of interest; once found, observatory slews towards source and uses second 
instrument with narrower field of view/higher resolution (e.g. Fermi, Swift) 


— (2) Surveys are planned prior to mission launch or during the mission to point at 
portions of the sky where known targets of interest reside (e.g. HST, JWST) 


— WEIRST falls into the latter category 


¢ Determining the worst thermal cases for survey type missions can be challenging given 
the large range of pointing possibilities 


— Worst case may be at edges of Field of Regard, but not necessarily 
— Determining worst case slew for stability requirements is also challenging 


— ls worst case slew from one Field of Regard orientation to another realistic, based on 
the expected mission profile, especially for known target surveys? 


¢ WFIRST has representative pointing orientations for entire mission to ensure that all 
Surveys can be completed in mission lifetime with allowances for guest observer time 


— These variations and dwells may occur over short time scales 
— Not every slew is thermally significant. How can key changes be identified? 
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The entire sequence of pointings for WFIRST was studied by the systems 
engineering team to seek a single 24 hour period which featured frequent and 
extreme slews as a “Day-in-the-Life” analysis case 


This resulted in 388 slews on approximately a 4 minute cadence. This was deemed 
to be too many points to recalculate the thermal environment each time 


Need to identify significant changes YPR Curve Fit for TD parameters 


Roll --¢-- Yaw Fit --¢-- Pitch Fit —+— Roll Fit 


Pitch 











Yaw 


— Initially tried to use angle between = 
sun vector for subsequent points (dot 
product). Did not account for sign... 








‘ ; ; ® 100 
— Used variation in vector component Bp 
C 80 
— Average X, Y, and Z component of S . 
normalized vector to sun since last - 
calculation point computed a eee a fit 
20 
> | 
— When any parameter deviated by a J, | 
more than 10% since last calculation 20 
point, significant change identified 0 
0 7200 14400 21600 28800 36000 43200 50400 57600 64800 72000 79200 86400 
— 36 Calculation points identified Time (s) , 


a ee an 
Ca i 


mVCHUr aaa aD, 


WIDE-FIELD INFRARED SURVEY TELESCOPE 
ASTROPHYSICS * DARKENERGY ¢ EXOPLANETS 





1 
%, ; 4 \ / ; . : \ yy i ; 
vA yy * te) oe 
F ne iy ea 
a AN og 
Te \ yd Dat). 
- If st ni 7 My if Vv 
We Way) HL 
Cae a ot \ 
; vr We ee PAA) WAR 
i ’ "yh Vy iv 4 ih \ § ta ai ra 
i” i A 4 
; Vola M ve hy a . e 
“§-" 4 iy. ah ah 
lly NANG as!) : y ‘ 
' NU ide 
: i  \oo 1 
"Ne L yy 
Ae 4 
1g ek MOND v 
eye At » WR 
‘ 4 A MG > } 
: , wR ie ‘Vy , 
, en ‘> 
\ = & y i 


¢ WFIRST has a number of PID or PI heater controllers to meet stability requirements 
¢ Predicts for one controller showed apparently poor control. 
—In actuality, the non-uniform output sampling frequency was misleading 


¢ Even with uniform output frequency, the stability requirements were still not met 





. Q - p * (T —T ) + 982 1 ou bSample:Sensing ---- FullSample:Sensing —*—SubSample:Heat  ---- Fullsample:Heat ' 
HTR GAIN SP ACH 


Iocan & (Tsp — Tacy) * at + _ 








DAN * (Tsp 7 Tacu)/ dt Temperature | 
¢ Initially believed that P.,,,, was VA i el & 


__ 282.02 : $s \ yi LN i . | y : : , : 0.6 
too high ¥ Rk AAAS MOY PE WN AVE RIT | 2 
2 22 YH RRS ese eee so) 
— Minimal impact on performance’ _ 2 z 
281.98 0.4 
when P-.,,,, was chan 
En Pgain Was changed — Heater Power 


¢ How can the gains be adjusted PA ga ga 
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° ° ° ontributions from P an omponents to lotal Heater Power 

¢ PID modeled with prevention of — aie tee <ceme Sones 

windup and constraints on the heater power _ ye 
— Long warmup or cooldown periods can =. — 
result in long term error sum that takes z “ou 
many cycles to eliminate 8 22.002 
— Heater had min of 0.0 and max of 0.25 W : _ a: 

¢ Control Temp output at each timestep 5 _— 

° Based on the Error and the Error Sum, the eae 76200 76400 76600 76800 77000 77200 77400 — 

. . Time for One Heater "On" Cycle (s) 
contributions from the PAIN and GAIN were Comparison of Original Response to Response with Reduced lean 
calculated individually _ —trgedanete. —iedgacintaseiee 

— These values were plotted over a time when 282.1 

the heater was active to identify which 282.08 
component was the driver (Top Right) F 282.06 

—_— 

— The le, clearly is shown to be the larger oem 
contributor aes 

¢ Adjusting the I.,,, resulted in the much _ 


better control (Bottom Right) 0 43200 86400 129600 172800 216000 


Time (s) 8 
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¢ Laptop used for runs: Intel® 7-37200M with the CPU running at 2.6 GHz 
—4 Cores, 4 Virtual Cores with Hyper Threading 


— 32 GB of RAM — no need for Virtual Memory j 5 A 8 
SINDA/FLUINT v5.7, Patch 9 Time to Run Single Job (hr) 

¢ Compare Field of Regard and Slew run times 942 9 39 1212 19.47 

—1, 2, 4, and 8 jobs submitted in parallel 17.32 17.85 20.48 31.85 


8 x Time(1 Job Serially) Time to Run 8 Jobs (hr) 
75.33 39.27 24.23 19.47 


138.53 71.40 40.97 31.85 


—Time(8 jobs serially) = 
e Exact same model so results are identical 


¢ Running 8 jobs in parallel requires about 
twice as much time to get any one set of Periormance Degradation 
results for Single Job 
. . 0% 4% 29% 107% 
¢ Running 8 jobs serially takes about 4 times 
. 0% 3% 18% 84% 
longer than 8 jobs in parallel to get all results 
A Jobs in Parallel | f Performance Degradation 8 Jobs 
° in Par isa “sweet spot” for gettin 
obs aralle isa S ee spo | or getting Seren ane (ae a 
some results earlier without waiting too long 
335% 124% 29% 0% 
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¢ Fluid modeling coupled with a thermal model adds complexity due to differences in 
allowable timesteps based on the scales of the thermal domain and fluid domains 


¢ Time dependence effects on the Fluid Lumps and Paths impacts run time (Junctions and 
STubes are time independent, Tanks and Tubes are time dependent) 


¢ Accuracy impact for this model is considered negligible 
¢ Surprisingly, Junctions (Time Independent) and Tubes (Time Dependent) resulted in fastest 
run times. 


— Hypothesize that Tubes allow slightly larger timesteps and may have reduced backup tries when 
convergence criteria not met 


Min Err Avg Max Err Avg Max Err 
(K) Err (K) (K) Min Err (K) Err (K) (K) 
10.92 0% 2535 25416 -- -- ~~ -- -- ~~ 
7.70 -29% 2242 21223 -0.486 0.002115 0.342 -0.197 0.004625 0.091 
24.002 120% 10102 #£=56031 -0.190 0.001614 #£1.057 -0.174 0.004452 0.077 


31.63 190% 14297 70866 -0.450 0.002363 0.990 -0.174 0.004663 0.077 
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¢ The overall model runtime is sively ifiuasua by the size of the solution matrix 
¢ Density of the solution matrix strongly related to the number of radiation couplings 


— Radiation coupling filtering can substantially reduce the density of the matrix, but this 
filtering is generally based only on the interchange factor and not the heat exchange 


¢ Time to fire more rays can also impact overall run time 
— Many thermal cases might use one set of radiation results 
—MCRT is inherently parallelizable 
¢ Regardless of Number of Rays fired, number of output radks is nearly the same 


— Not worth the time to fire more rays and then throw away most of the small radks 


0 0.001 0.9 14.26 6.55 20.81 6,330,895 16,595,983 
1 0.001 0.9 14.25 6.4 20.65 6,332,632 16,601,196 
0 0.001 0.9 Ios 6.07 66 6,257,883 31,678,395 
1 0.001 0.9 60.6 6 66.6 6,262,274 31,671,957 
0 0.001 0.9 608.4 9 617.4 6,137,361 58,939,273 
1 0.001 0.9 575.4 7.8 583.2 6,140,926 598,396,535 
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Filtering used to output more 
or fewer radiation couplings 
from the solution database 


Filtering is performed by 
sorting all Bij terms from high 
to low and including all terms 
above Cutoff 


Including a Sum term will 
add in additional Bij terms 
until the Sum criteria is met 


Again, regardless of the 
number of rays fired, the 
total number of radks output 
is about the same for all 
other parameters being equal 
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0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 


0.95 
0.95 
0.95 
0.95 
0.95 
0.95 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
N/A 
N/A 
N/A 
N/A 
N/A 
N/A 


10,142,081 
10,152,226 
10,290,886 
10,285,600 
9,705,460 
9,702,782 
6,330,895 
6,332,632 
6,257,883 
6,262,274 
6,137,361 
6,140,926 
2,748,845 
2,/48,614 
2,681,166 
2,680,655 
2,668,267 
2,668,210 


12,760,169 
12,757,101 
27,622,053 
27,625,250 
54,951,103 
54,814,741 
16,595,983 
16,601,196 
31,678,395 
31,671,957 
58,539,273 
58,396,535 
20,199,716 
20,207,021 
35,273,939 
35,272,338 
62,025,499 
61,886,396 


22,902,250 
22,909,327 
37,912,939 
37,910,850 
64,656,563 
64,517,523 
22,926,878 
22,933,828 
37,936,278 
37,934,231 
64,676,634 
64,537,461 
22,948,561 
22,955,635 
37,955,105 
37,952,993 
64,693,766 
64,554,606 
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¢ WFI Instrument includes a re wdiels 


provides some key metrics for comparing N/A 0.9 0.95 
predicts 5.78 12.7 20.73 
¢ Three runs compared: 
6.712 8.035 9.442 
— Cutoff = 0.001, No Sum, 2.7 M Radks, 5.8 hrs 945 112.0 145.2 
— Cutoff = 0.001, 0.90 Sum, 6.3 M Radks, 12.7 hrs 734.8 J48.8 971.5 
— Cutoff = 0.001, 0.95 Sum, 10.1 M Radks, 20.7 hrs Avg Err AvgErr AvgErr 
¢ Model with 0.95 Sum took nearly 3.5x longer (K) (K) (K) 
than No Sum, but accuracy is questionable... — ie eas 
’ -- 14.19 22.88 
¢ With 0.95 sum, load on CC nearly 40% higher!! = 3244 5.03 
— This does not make sense... a 0.06 -1.26 
. ; - -- -1.84 -3.64 
¢ Impact is felt on 270 K Radiator, which is 
rejecting Cryocooler Compressor power = bated a 
JECHINE MY p p z 1.21 2.28 
¢ IFU Detector is also greatly influenced a 7.94 19.39 
-- -0.06 -0.14 


¢ Grism and Filters are also influenced, but effect 
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Why is CC Load so different? Can 5-10% 
more FOV really cause this? 


should be small on these components 
¢ Further investigation needed... 
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118.5 0 110.55 0 99.112 0 
1.576 1.614 1.647 
-2.667 -2.321 -1.938 
1.081 0.7096 0.2853 
0 0 0 
Type Tempj Cond} Heat | Temp j Cond j Heat | Temp j Cond j Heat | 
L 92.73 21.395 -2.667 92.589 21.395 -2.321 92.419 21.395 -1.938 
L 165.55 2.52 1.576 165.75 2.52 1.614 165.87 2.52 1.647 
R 163.25 5.15E-07 1.80E-05 164.6 5.15E-O7 1.88E-05 165.94 5.15E-07 1.96E-05 
R 164.78 1.71E-05 0.0001 165.03 1.71E-05 0.0001 165.15 1.71E-05 0.0001 
R 165.99 0.0001 0.0008 166.18 0.0001 0.0008 166.3 0.0001 0.0008 
R 164.41 0.0009 0.0314 164.94 0.0009 0.0322 165.39 0.0009 0.0328 
R 166.75 0.0024 0.0657 166.69 0.0137 0.4286 0 0 0 
R 165.55 0.003 0.0718 165.75 0.003 0.073 165.87 0.003 0.0741 
R 166.06 0.0057 0.1724 166.24 0.0057 0.175 166.37 0.0057 0.1775 
R 164.63 0.0288 0.739 0 0 0 0 0 0 


¢ Grism and Filters should not even be able to “see” the IFU Detector... 
¢ Why does adding in more couplings drastically change the heat flow to Grism and Filters? 
¢ Same effect seen for IR MOSAIC plate...software vendor contacted 
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¢ Software vendor found an error “7 


the Summation routine related to N/A 0.9 0.95  N/A::0.95 
5.87 12.32 18.18 10.38 
the model size and addressable 
memory 6.72 6.91 7.10 7.09 
. 9459 96.05 97.29 97.17 
. Quickly fixed and results regenerated 234.2 236.1 237.3 237.2 
¢ Errors were much more acceptable ne oa — 
¢ Run time with 0.95 Sum was about 0.0286 0.0113 : 0.0316 
0.8293 0.2939 = -0.6533 
50% more than 0.90 Sum and 300% Diss cise Bees 
more than No Sum 0.1034  -0.0067 = 0.1689 
— 0.0452 -0.0384 = 0.1015 
e Errors were generally within 1K 0.3786 0.1736 i 0.0087 
0.3635 0.1634 = 0.0081 
e@ 
Hybrid approach developed which ne re = caaES 
extracted Radks from 95% sum case -0.0490 -0.0278 Es -0.0143 
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for cryo components and No Sum 
case for non-Cryo components 


¢ Cryocooler is sensitive to more of the Radks 
that are filtered than other regions 

— This results in better prediction of the 
Cryocooler load while still improving 
run time over 0.90 Sum case 


¢ Would be a nice feature for S/W vendors to 
filter based on Energy (Bij * (T* — T*)) 
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¢ WFIRST is a large and complex model already, even in Phase A 
— Complexity driven by requirements and analysis questions (STOP, tight stability) 
¢ As Phase A Trade studies are underway, analytical efforts are needed to judge and evaluate 
potential designs and to verify that requirements can be met 
¢ The ability to quickly exercise the model is critical to providing data 
— Simulating realistic slew profiles instead of only Field of Regard constraints 
— Ability to reasonably adjust and tune PID parameters to meet stability requirements 
¢ Effective usage of computational resources 
— Parallel job submission should take into account if partial data is needed sooner 


— Fluid modeling with time dependent lumps greatly increased run time, but seemed to decrease run 
time for time dependent flow paths 


— Firing more rays and filtering the small terms wastes computing time 
— Including more (smaller) radks increases run time, but accuracy may be necessary for cryo regions 
¢ Filter Radks based on Energy from previous runs results instead of just geometric factors 
¢ No substitute for careful evaluation of model predicts (models must still follow physics!) 
— Bug found in code that gave questionable results. Correction resulted in more realistic predicts 
— The software will not tell you when it is wrong and cannot replace analyst’s judgment and experience 


¢ WFIRST now in Phase A and proceeding to SRR...more to come next year! 


